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ABSTRACT 

Aims. We improve the description of the evolution of the Sun's open and total magnetic flux on time scales of years to millenia. 
Methods. In the model employed here the evolution of the solar total and open magnetic flux is computed from the flux emerging 
at the solar surface in the form of bipolar magnetic features, which is related to the sunspot number cycle parameters and can be 
estimated from historical records. Compared to earlier versions of the model in addition to the long-lived open flux, now also a more 
rapidly decaying component of the open flux is considered. The model parameters are constrained by comparing its output with 
observations of the total surface magnetic flux and with a reconstruction of the open magnetic flux based on the geomagnetic indexes. 
A method to compute the Sun's total magnetic flux and the sunspot number during the Holocene, starting from the open flux obtained 
from cosmogenic isotopes records, is also presented. 

Results. By considering separately a rapdly evolving and a slowly evolving component of the open flux the model reproduces the Sun's 
open flux, as reconstructed based on the aa-index, much better and a reasonable description of the radial component of interplanetary 
magnetic field data are obtained. The greatest improvement is in the reproduction of the cyclic variation of the open flux, including 
the amplitudes of individual cycles. Furthermore, we found that approximately 25% of the modeled open flux values since the end 
of the Maunder Minimum are lower than the averaged value over 2008, i.e. during the current low minimum. The same proportion is 
observed in reconstructions of the open flux during the Holocene based on cosmogenic isotopes, which suggests that the present solar 
minimum conditions are below average, but not exceptional in terms of the heliospheric magnetic flux. 

Key words. Sun: magnetic fields, Sun: sunspots, Sun: solar wind, activity, Sun: evolution, Sun: solar-terrestrial relations 



1. Introduction 

The variability of the magnetic field has a strong influence on the 
dynamics of the outer layers of the Sun. Thus, the 1 1 -year cyclic 
variability of the magnetic field is registered by several solar pa- 
rameters such as the sunspot number and area, the rate at which 
flares and coronal mass ejections occur, the flux of solar X-rays, 
radio waves and solar energetic particles, as well as the total and 
spectral solar irradiance. The variation of the total magnetic flux 
and the surface distribution of the field also influences the open 
magnetic flux and hence the heliospheric magnetic field. 

Our knowledge of the evolution of the Sun's magnetic flux 
on longer time scales is limited by the availability of continu- 
ous and reliable observations of the solar magnetic field. Most 
of these observations are available for just a few decades, e.g. 
since the beginning of the space age in the case of the open mag- 
netic flux. On longer time scales, the solar magnetic flux must be 
reconstructed or computed from proxies. Thus, the heliospheric 
flux (i.e. the solar open flux) is reconstructed b ased on the ge- 
omagnetic aa-index from 1 868 to the present dLockwood et alj 
119991: iRouillard et alJl2007» . The open and total magnetic flux 
since roughly 1610 c an be computed from the sunspot number 
dSolankietalJl2000l |2002|) . During this period, the open flux 
doubled. On millennial time scales, reconstructions of the so- 
lar open flux and sunspot number were obt ained based on cos- 
mogenic isotopes such as 14 C and 10 Be ( Sol anki et alj l2004t 
Usos kin et al"1l2002Ll2003l. 12004 120071) . 
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The evaluation of the solar surface magnetic flux bud- 
get requires continuous and full-disk magnetograms, which 
are available for just for the last few solar cycles. Even 
for this period the total amount of magnetic flux emerg- 
ing in small bipolar magnetic regions is uncertain due to 
the li mited spatial resolution of the data dKrivova & Solanki 
|2004|) . In order to assess the solar magnetic flux budget on 
longer time scales s emi-empirical models hav e been used 
dKrivova et alj 120071: ISolanki et al] 120001 |2002|). as well as 
flux transport computatio ns (Baumann et al. 2004; M ackav et al.l 
l2002l:IScnussler & Baumannl2006HWang et al.l2002ll2005h . The 
former are based on the sunspot record and attempt to recon- 
struct the evolution of the total magnetic flux and the flux emerg- 
ing in large and small bipolar regions as well as the evolution of 
large unipolar regions that give rise to the solar open flux. These 
models are validated by comparing them with the record of the 
total magnetic flux obtained from magnetograms as well as with 
a longer record of the Sun's open magnetic flux deduced from 
the aa index, a measure of the variability of the Earth's magnetic 
field produced by its interacti on with the variable in terplanetary 
field, i.e. the Sun's open flux dLockwood et alJfl999h . 

The models of lSolanki etail d2000l[2002h of the evolution of 
the solar magnetic flux reproduce the long-term variation of the 
solar open flux reasonably well, but give a much smoother vari- 
ation over the solar cycle than the observations. In addition, the 
computed open flux lags the observations by roughly 2-3 years. 
These shortcomings are the result of the simplicity of the model, 
which, in order to explain the long-term trend requires the open 
flux to have a lifetime of multiple years. Such an extended life- 
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Fig. 1. (a) Comparison of the observed (symbols) and modeled (solid line) total magnetic flux. Each data point is an integral over 
a synoptic chart for one Carrington rotation. Different symbols are used for different data sets: circles represent the KP NSO data, 
squares MWO data and diamonds WSO data. For the modeled flux the value (<p ac , + 03<p ep k + (p D pen) is given (see text). The dashed 
blue lines bound the period used for the optimization of the parameters of the model between 1974 and 2002. (b) Reconstructed 
magnetic flux in AR (blue line), ER (green line), open (Red Line) and total flux (cyan line). For the total flux here the value 
((pact + <f>eph + ^open) is plotted, (c) Fractional contribution of AR, ER and open flux to the total Flux. 



time is appropriate for the field in the large polar coronal holes, 
but not for some of the smaller, relatively short-lived low-latitude 
coronal holes that are the source regions of a significant fractions 
of the open flux during the high activity phases of solar cycle. 

The main motivation of this work is to estimate the evolution 
of the solar magnetic flux taking into account that in addition 
to the long-lived open flux there are also more rapidly evolving 
source regions (coronal holes). These are often associated with 
active regions or decaying active regions. For the purposes of 
the model, it is important to distinguish between the slowly and 
rapidly evolving flux and not between the locations on the solar 



surface. In particular, there is no one-to-one correspondence of 
high and low latitude coronal holes to slowly and rapidly evolv- 
ing flux. 



The paper is structured as follows. In Sect. 2 we describe 
the solar magnetic flux model, the parameter optimization pro- 
cedure, and the model parameters. The results are described in 
Sect. 3. The extension of the model to periods prior to telescopic 
sunspot observations is presented in Sect. 4. The conclusions are 
given in Sect. 5. 
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2. Model Description 

2. 1 . Approach 



bv lSolanki et al.l (I2002I) . which is itself an extension of the work 
of ISolanki et alJ (l2000h . In this model the evolution of the to- 
tal and open magnetic flux is computed from the flux emerg- 
ing at the solar surface in form of bipolar magnetic features. 
On the Sun, th e bipolar mag netic regions display a continuous 
size spectrum (iHarvevlll993l) . In the model the spectrum is di- 
vided into two classe s accord i ng to t he size and life-time of the 
structures, following lHarveyl (Il993l) and others. Large bipolar 
structures emerging in the activity belts and living up to several 
weeks are classified as Active Regions (AR). Small short-lived 
bipoles, which emerge over a larger range of latitudes, are clas- 
sified as Ephemeral Regions (ER). The magnetic flux emergence 
rate in active regions is roughly proportional to the sunspot num- 
ber, which allows it to be estimated from historical records. The 
emergence rate of magnetic flux in ephemeral regions is higher 
than in active regions and their contribution to the total photo- 
spheric magnetic flux is significant. In spite of their limited life- 
time the number and latitude of ephemeral regions also evolves 
over a cycle, which is extended with respect to the sunspot cycle, 
but shows a much smaller contrast between activity maximum 
and minimum. The literature contains contradictory statements 
wheth er there are more or less epheme ral regions at activity max- 
imum (Hagenaa3|2001tlH arveyll 1993b . In the model, parameters 
of the cycle (such as time of maximum, amplitude and length) 
displayed by the ephemeral regions are assumed to be related to 
the properties of the corresponding sunspot cycle. 

Part of the magnetic flux that emerges in active and 
ephemeral regions is dragged outward by the solar wind and 
reaches far into the heliosphere. It is called the open magnetic 
flux. As the source of the open magnetic flux is located in of- 
ten large regions with a dominant magnetic polarity, it can sur- 
vive on the solar surface for a relatively long time, reaching 
up to several years. There are, however, also smaller, shorter 
lived coronal holes often associated with decaying active re- 
gions. These lead to a far more rapid variation in the Sun's 
open magnetic flux, in particular around activity maximum and 
shortly after it. Hence, some of the flux from active reg ions does 
open, but stays open only for a r elatively short time (ICranmerl 
2002). [ikh sanov &Ivanovl (119991) show a histogram of equato- 
rial coronal holes (CH) lifetimes according to which nearly 50% 
do not outlive 3 solar rotations, implying a median lifetime of 
80-90 days. This may be an upper limit, since they only consider 
equatorial coronal holes that survive at least 2 solar rotations. We 
stress, however, that our rapidly and slowly decaying open fluxes 
cannot be simply associated with low and high latitude coronal 
holes. An alternative interpretation of the rapidly evolving open 
flux is that it is additional flux carried into the heliosphere by 
CMEs before their disconnection from the Su n ( Crooker et al.l 
l2002t iLuhmann et allll998t lOwens & Crookerl 12006). As such, 
it is not strictly open (in the sense that it does not reach all the 
way out to the Heliopause), but it does contribute to the inter- 
planetary field and the Sun's magnetic flux at 1 AU. This is a 
relevant quantity for comparing with th e open magnetic flux at 1 
AU reconstructed by iLockwoodl d2009l) . 

In the followin g, we extend the previous model of 
ISolanki et al.l (|2002) by distinguishing between rapidly decay- 
ing open flux, <p r open , and slowly decaying open flux, <p s open . A set 
of four coupled ordinary differential equations then describes the 



evolution of the four surface magnetic flux components that we 
consider: 



d<p a 



Th e model presented h ere is an extension of the one presented dt 



Fad 

To - 



"act 



fact 

~ 

' ai l 



d<peph 

dt 



<Peph <Peph 
Ceph ~ "J -j— 



d<p, 



open 



dt 



dt 

bopen 
btotal '■ 



fact 



act 



eph 

$apen 

T r 

open 

4>eph 
eph 



eph 



1 open 



Popen $open ■ 
f>act + <Peph + < 



'open • 



(i) 

(2) 
(3) 

(4) 

(5) 
(6) 



where 4> act , <p ep i t , cpopen, and 4>totai refer to magnetic flux of AR, 
ER, open flux, and the total flux. The open flux is the sum of the 
rapidly (<p open ) and slowly (<p s open ) evolving components (Eq. [5}. 
Distinguishing between these two components of the open flux is 
the ma in difference to the c omp utations of magnetic fl ux evolu- 
tion bv lSolanki eTai] d2002l) and lKrivova et alJd2007l) and leads 
to the introduction of one more differential equation. The <p r open 
is fed only by <p act since it is assumed to reside in small coronal 
holes located close to ARs. In contrast, <p s open obtains contribu- 
tions from both active regions and ephemeral regions. It resides 
at least partly in the polar cap coronal holes, but also elsewhere 
(e.g. in the flux from decaying active regions wandering to the 
poles). 

The time constants t° c/ , T° ph , T r open and r s opm are the decay 
time scales of the AR, ER, rapid and slow compone nts of the 
open flux, respectively. Following Solanki et al. (2000) and oth- 
ers, we assume that the decay process is due to the cancellation 
with flux of opposite polarity, but do not specify the process be- 
yond giving the decay time. The time constant T 2 act is the flux 
transfer time from active regions to the rapidly evolving com- 
ponent of the open flux, while T l act and r l eph are the flux transfer 
times from AR and ER to the slowly evolving component of the 
open flux, respectively. 

The input parameters of the model are the flux emer- 
gence rates of active and ephemeral regions. Here, the flux 
emergence rate in active regions is chosen to be linearly 
proportional to the mo nthly averaged group sunspot number 
dHovt & Schattenll 19981). Ra, and is sca led according to the ob- 
servations of[S chriiver & Harvevld 19941) for cycle 21. Following 
iKrivova et al. (2007), we define the flux emergence rate of AR 
as 
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with e™' v ' 21 = 2.3 x 10 24 Mx yr" 1 andfl"""' 21 = 172. The value 



of R™ ax < 2x was obtained from 3-mon th running means of R g . 

Following ISolanki et al.l d2002l) . we obtain the total flux 
emergence rate in ER as a sum over multiple overlapping cycles 
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where we define e' eph (t) as 

i Ph {t)=crxg\t). (9) 

Here X is a scaling factor and g'(t) is a function defined as 



As the error of individual observations is unknown, we use 
the standard deviation of the observed <p to tai and reconstructed 
<t>one. n data sets to estimate th e values of x' 2 ■ Following oth- 
ers dHolland & Welsch"lll977l) . we apply a weighting function 
(Cauchy weighting function) 



g\t) = 



COS' 

,0, 



-T' eph /2<(t 
other wise. 
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(10) 



1 +r 2 
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where t' c is the time at which cycle ; reaches maximum activity 
and T' h is the length of the ephemeral cycle. Equation (TlOb is 
equivalent to Eq. (6) of iRrivova et al. (2007). The length of the 
ephemeral cycle is related to the length of the activity cycle, T', 
in the following way: 



T' -T l - r 

1 ext 1 c -» ' 

and 

7^ = ^+27^, 



where c x is the ER cycle extension parameter and 2T' exl is the 
extension of the ER cycle i in relation to AR cycle i. 

We define the maximum emergence rate at cycle i, which 
appears in Equation (O as 



jnaxd max,2l 



r: 
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where R g 
i. 



is the maximum sunspot number observed at cycle 



2.2. Parameter optimization 

The model has 8 free parameters. These are introduced and dis- 
cussed in Sect. 2.3. Here we describe how the optimization is 
carried out. The free parameters were adjusted by comparing 
the model output with observations of the total surface mag- 
netic flux deduced from synoptic charts of the Sun's radial field 
( Arge et alj |2002|) and with a reconstruction of the open mag- 
netic flux based on the geomagnetic aa-index by iLockwoodl 
(120091). This is a revised version of the reconstruction due to 



Lock wood et alj (1 19991) . This revised open flux includes correc- 
tions due to kinematic effects produced by the propagation of 
CMEs (Lockwood et al. 2009a, b). It covers the period from 
1904 to 2008. 

The optimization of the model's free paramet ers is realized 
using the genetic algorithm PIKAIA described by Charbo nneaul 
( 1995). The code maximizes a function (fitness function, /), de- 
fined as 



f {parameters) 
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X 1 total ^ X' open 



(14) 



where x' 2 is the reduced^/ 2 , i.e. the^ 2 per degree of freedom 
(Df). For this analysis, we define 



*' 2 = ^2> 

J 1=1 



obs 



In order to reduce the influence of outliers, the value r in the 
weight function is set to 



'( 



ad(h) \ 



(17) 



0.6745 ; 



(11) where 9 is a tuning parameter, h is the deviation of the model 
' , ' ). Here 

■s' 



model 



from the observations/reconstruction (h 

mad is the median absolute deviation of the residuals from their 
median values, and the constant 0.6745 makes the estimate un- 

(12) 

biased for the normal di stribution. The tuning p arameter (0) is 
set to 2.385 according to lHolland & Welschl d 19771) . 



where /V is the number of observed data points, while x 
and x' , , represent the i-th observed and modeled data point, 

model r r 

respectively, cr' is the error of the i-th observed data point, tion of < 



2.3. Input data and parameters of the model 

For the total surface magnetic flux we use a set of ob servations 
compiled bv lArge et all (120021) and lWang etaTI d2006l) . This data 
(13) set is based on almost daily observations of the solar global pho- 
tospheric field that have been carried out at the Mt. Wilson Solar 
Observatory (MWO), National Solar Observatory Kitt Peak (KP 
NSO), and Wilcox Solar Observatory (WSO) over cycles 20-23. 
Note that for the optimization we only use data recorded be- 
tween 1974 and 2001 when data from all observatories is avail- 
able. Outside these periods we compare with the observations a 
posteriori as a further test of the model . We take into account 
the finding of iKrivova & Solankil d2004l) that more than half of 
the photospheric flux of ER may escape detection in the em- 
ployed synoptic maps due to their relatively low spatial resolu- 
tion. Consequently, we compare the measurements of the total 
magnetic flux to the value (<jv, + c eph $ eph + (popen), where c eph 
is the fraction of <p ep i, that is detected, which is relatively un- 
certain, because magnetic polarities are often missed on small 
scales. Here, we set c epn to be approximately 30% c ompared to 
the value of 40% employed bv lKrivova et alj (120071) . The value 
of c ep h affects the amplitude of the ER cycle in the model since 
the values of the total flux during solar activity minima is deter- 
mined mostly by the ER flux. 

The parameters of the model as well as their adopted or best 
fit values are listed in TableQ] Th e ER flux decay time is fixed to 
a valued of 14 hours as found bv lHagenaarl d2001l) . The remain- 
ing parameters are allowed to vary within a given range based 
also on independent observations and/or physical assumptions. 
This range is also given in Table Q] 

The ranges of the decay time scale of the AR flux, r° ct , 
the ER amplitude factor, X, and th e extension pa r amete r, c x , 
are set as di s cussed prev iously by IKrivova et alj ([2007) and 
ISolanki et all (l2000i |2002|) . The best value found for the decay 
time scale of the AR flux is approximately 0.32 years. This es- 
timate is close to the 0.25 years obt ained assuming a b alance 
(15) between flux emergence and decay (IKrivova et al.ll2007h . The 
best-fit ER amplitude factor is approximately 106 and the cycle 
extension parameter is approximately 5 years, which leads ER 
cycles that are longer than 20 years. These values are different 
from previouly obtained ones, which indicates that the introduc- 



changes the best-fit solution. 
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ISolanMetaTld2002l) and lKrivova et al.l(l2007h postulates that 
the flux transfer time scale from ER to the open flux is a fac- 
tor of six higher than the flux transfer time scale from AR 
(j eph = 6T l ac[ ). This assumption is based on the observation that 
in cycle 21 the contribution of ER to the axial dipole moment of 
the Sun was about a factor of six smaller than t hat of the AR , 
assuming an average lifetime of ER of 8 hours (lHarvevl [T994). 
Here, we have not constrained the value of T eph in relation to 

T act because it is not clear if only T s open contributes to the ax- 
ial dipole or if T r open also influences it. In addition, we employ 
a longer average life time of ER (14 h ours) as proposed by the 
more recent work of lHagenaarl d200ll) . We searched the best-fit 
values of both parameters in a wide range from 10 to 90 years. 
We defin ed this wide range in order search values in the domain 
in which Kri vova et alj d2007l) found the solutions for this pa- 
rameters. The best fit values for T l act and t h are approximately 
85 and 10 years, respectively. 

Note that due to the extended length of the ephemeral cycle, 
around activity minimum both the preceding and following cycle 
contribute to the total and open flux. For the current minimum, 
we do not yet know the features (strength, length, and time of 
maximum) of the next cycle (cycle 24). Therefore we negled 
this cycle completely, so that the modeled magnetic flux values 
during the current minimum may be too low. 

We searched for the transfer time of AR flux to the rapidly 
decaying open flux (T 2 act ) in the range between T° eph (14 hours) 
and 6 years. The best-fit value found for this parameter is ap- 
proximately 1 .7 years, while the best value of the decay time of 
the rapidly evolving component of the open flux (T r open ) is about 
46 days. 

The decay time scale of the slowly evolving component of 
the open flux (j s open ) is restricted to the range between 50 days 
and 6 years. The upper limit was defined in order to contain the 
values previously found for the decay time scale of the open flux 
(approximately 3-4 years). The optimum value returned by the 
code is approximately 1.4 years. This value is shorter than the 
previous estimate of the decay of the open flux between 3 and 4 
years. 

We note that the set of best-fit parameters presented in Table 
1 do not constitute a unique solution. In particular, we found so- 
lutions for short values of T l eph ( 3 years) and T s open ( 7 months) 
that have similar values of the fitness function. While a de- 
cay time of t he open flux of about 3-4 year s, as found by 
IKrivova et all d2007l) : ISoranki et all d2000ll2002l) . determines the 
long-term evolution of the open flux by producing significant 
overlap between cycles, a short decay time of approximately 7 
months cannot on its own account for the observed secular vari- 
ations because the flux from the previus cycle has decayed be- 
fore the next cycle starts properly. Instead, for such a short t s n 



most of the long-term variability of the slowly evolving compo- 
nent of the open flux is due to the flux transferred from ER. For 
T open - 1-4 years both process play a role. 

3. Results 

3. 1 . Comparison of the model output and the 
observations/reconstructions 

3.1 .1 . Total Magnetic Flux 

Figure [T^ shows a comparison of the observed (symbols) with 
the modeled (solid line) total magnetic flux. The three obser- 
vational data sets are represented using different symbols: cir- 
cles represent the KP NSO data, squares MWO data, and di- 
amonds WSO data. Values are given for each Carrington ro- 
tation (CR) from the start of observations to the present: CR 
1615-1975 (NSO), CR 1516-2082 (MWO), and CR 1642-2081 
(WSO). Only the period between the vertical dashed lines is used 
for the optimization of the parameters. The total flux plotted here 
is given by (p, ota i = <p act + 0.3(p eph + 4> open . The factor 0.3 takes 
into account that a major part of the photospheric flux from ER 
is missed when emplo ying synoptic charts due t o their relatively 
low spatial resolution dKrivova & S olanki 20041) . 

The model reproduces well the average variability of the data 
sets employed. The minima between the cycles 21-22 and 22-23 
are slightly overestimated while the minima between the cycles 
20-21 and 23-24 are well reproduced. Note that the observations 
of the descending phase of cycle 23 and the minimum between 
the cycles 23-24 were not employed for the optimization of the 
parameters of the model. It is gratifying to see that the model 
reproduces the total magnetic flux during the cycle 20 and during 
the declining phase of cycle 23, although these data points were 
not used to constrain the solution and the low value of the flux 
during the current minimum lies well outside the range of the 
previous minima. 

Figure lb displays the calculated evolution of the AR mag- 
netic flux, ER flux, open flux and the total flux. In this panel the 
total flux is given by the expression: <p, ota [ = (pact + <peph + 4>open- 
The modeled minimum values of the total flux between the cy- 
cles 20-21 is reduced by approximately 20% relative to the min- 
imum value between cycles 21-22 and it has the same level 
as the minimum between the cycles 22-23. The modeled value 
of the present minimum is approximately one third of that fol- 
lowing cycle 22. For the model discussed here the ER flux has 
lower variation over the a ctivity cycle than prev iously found by 
IKrivova et all d2007l) and ISolanki et all d2002l) . This is due to 
the very extended length of the ephemeral regions cycle, lead- 
ing to a larger overlap between them. This small variation of the 
ephemeral reg ion flux is in between the variations proposed by 
lHarvevI d 1993b and Hagenaarj d2001l) . Shorter ephemeral region 
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cycles, wit h less overlap , produce results t hat are more similar 
to t hose of iHarveyl (119931) . as modelled bv lSolanki etail d2002l) 
and lKrivova et al.l ([20075. Even longer ephemeral region cycles, 
with even more o verlap, produce maxima in (f> act in agreement 
with the results of Hagenaar (2001). 

The relative contribution of AR, ER and open flux to the 
total flux is presented in Figure lc. The average contribution of 
the AR flux to the total flux is approximately 50% during the 
activity maxima while the contribution from ephemeral regions 
is about 40%. During the minima, the ER flux contribution is 
about 80% while the contribution from AR flux is approximately 
10%. The open flux represents a small fraction of the total flux 
(approximately 10% during all phases of the cycle). 

3.1.2. Open Flux 

Figure [2^ s hows a comparis on between the open flux recon- 
structed by ILockwoodl d2009l) (blackline) and the modeled 
open flux (blue line) from 1904 to 2008. The square indicates 
the value observed in 2008 (average over the year), the tri- 
angle the modeled value of the open flux (for the same pe- 
riod). For c omparison, the green line represents the model de- 
scribed by iRrivova et al.l (2007). As the parameters of that 
m odel were obtained by f itting to the open flux reconstructed 
by lLockwood et af] d!999l) . we optimized the parameters of the 
model accordi ng to the procedu re described in Sect. 2.3, i.e. 
to the data of ILockwoodl d2009t) . prior to plotting its output in 
Fig. [2^. The parameters obtained for the iKrivova et al.l d2007l) 
model are presented in Table |2j where the limits within which 
the parameters were searched for are also given. A lag of ap- 
proximatel y 2-3 y ears between the open flux reconstru c ted by 
ILockwoodl d2009l) and that based on the IKrivova et aT] d2007l) 
model is evident. Phase shifts are also observed in flux transfer 
models d ue to the effect of the decay term in the flux tran sport 
equation dMackav et alj2002tlSchussler & Baumannl20 06). The 
calculated open flux evolution based on the present model repro- 
duces well the reconstruction of the open flux based on the ge- 
omagnetic aa-index. In particular, the offset betwe en model and 
observ ations, which is well seen i n the results o f Solanki et al.l 
(2002) and in the model based on IKrivova et al.1 (2007) is now 
gone. In addition, amplitudes of individual cycles are now in 
general better reproduced. For reference, a dashed red line rep- 
resenting the average value observed in 2008 is drawn. Clearly, 
the present model reproduces this valu e, while the previous , sim- 
pler version of the model employed bv lKrivova et al.l (|2007) does 
not. 

We note a discrepancy between the observed and modeled 
values for cycle 19, which is the strongest sunspot cycle ob- 
served since the Maunder Minimum. The cause of this discrep- 
ancy during the maximum of cycle 19 is not clear. It may be that 
for such a strong cycle, parameter val ues are diffe r ent f rom those 
valid for other cycles. For example. IWang et all (120051) require 
different meridional flow speeds from cycle to cycle in order to 
reproduce the Sun's open magnetic flux and to ob tain the polarity 
flip a t the poles from one minimum to the next (Bau mann et alj 
120041) . The meridional flow speed could influence our parame- 
ters, which we have maintained unchanged for all cycles. Also, 
other s olar parameters, suc h as the emergence latitudes, could 
differs (Solanki et al. 2008), which is not taken into account by 
our model. Alternatively, the reconstruction based on the aa- 
index could underestimate the peak value of the open flux dur- 
ing high solar activity. The sensitivity of high latitude stations 
to the auroral electrojet can be reduce d at high activity becaus e 
the electrojet drifts to lower latitudes ( Lockwood et al.l | 20 09c). 



This effect is observed, for example, in the apparent saturation 
or even decrease of the geomagnetic AE-index during intense 
magn etic storms ( Akasofu 1981; Feldstein 1992; Gonzale z et al.l 
1994). As discussed by lLockwood et all d2009d) . in principle, 
this non-linear effect is not significant for the aa-index because, 
although it is derived as a proxy of substorm activity, it is ob- 
tained from mid-latitude stations. Except during very intense 
magnetic activity, this effect is avoided for mid-latitude stations 
because the auroral elect rojet always migrates tow ard the station 
with increasing activity dLockwood et al. 2009c). We speculate 
that due to the exceptionally high activity during cycle 19 the aa- 
index may have underestimate the level of maximum magnetic 
activity all the s ame. 

According to Wang et al. (2006), the low-latitude component 
of the open flux closely tracks the Sun's equatorial dipole com- 
ponent, whose strength depends on the amount of flux present in 
the active regions and on the longitudinal distribution of the ac- 
tivity. Thus, asymmetries in the longitudinal distribution of large 
active regions can lead to large-amplitude variations of the low- 
latitude component of the open flux. As the model presented here 
does not describe the longitudinal distribution of the activity, the 
dips (Gnevyshev gap) observed in the open flux near the cycle 
maxima are not reproduced. 

As pointed out in Sect. 2, in principle the magnetic flux 
emerging in ERs belonging to cycle 24 contributes to the slow 
component of the open flux during the descending phase of cy- 
cle 23 and the present minimum. However, the flux emerging in 
ER for cycle 24 was not included since we do not know the rel- 
evant cycle parameters yet. Nonetheless, the model reproduces 
the present low level of the open flux. 

This result suggests that cycle 24 will be rather weak or peak 
very late. However, since we cannot be sure that the deduced set 
of parameters represents a unique solution, we hesitate to use 
the present model to quantitatively predict the strength or length 
of cycle 24. Test calculations indicate that even for a given set 
value of the model's free parameters only a function f(s, I) of 
the strength, s, and length, I, of the next cycle can be determined 
on the basis of our model. 

The evolution of the rapidly and slowly evolving components 
of the open flux is presented in Fig. 2b as well as the complete 
modeled open flux. Most of the cyclic variation is determined by 
the rapidly evolving component, while the slow component pro- 
duces a background field that varies from cycle to cycle, but only 
relatively weakly over a cycle. The slow component peaks gen- 
erally after the fast one in the descending phase of the sunspot 
cycle. The relative contributions of the rapidly and slowly evolv- 
ing components to the complete open flux are shown in Figure 
2c. A drop of similar magnitude as the current one (but starting 
from a higher level, so that the slow open flux did not reach such 
low levels as in 2008-2009) is seen between cycles 19 and 20, 
when a weak cycle followed a very strong one. During the so- 
lar minima, the open flux is maintained almost exclusively by 
the slow component (Figure 2c). According to the model, the 
rapidly evolving component contributes with about 40% of the 
open flux at solar maxima and this fraction remained almost con- 
stant during the last century, even during cycle 19. 

3.1 .3. Estimates of errors in the model 

Figure [3] presents further quantitative ways of comparing the 
model with the observations. Figure 3a displays the scatter plot 
of the modeled versus the observed total magnetic flux. As in 
Figure la, different data sets are indicated by different symbols. 
Each point represents a Carrington rotation. Also shown are a 
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Parameter 


Symbol 


Value (years) 


Min 


Max 


AR Flux decay time scale 


fact 


0.22 


0.2 


0.8 


AR Flux to Open Flux transfer time scale 


Tta 


18.46 


10.0 


90.0 


ER Flux decay time scale 


Teph 


0.0016 


Fixed 




ER Flux to Open Flux transfer time scale 


T,e 


110.76 


6r M 




Open Flux decay time scale 


t open 


4.87 


0.0016 


6.0016 


ER amplitude factor 


X 


159.05 


80.0 


160.0 


ER cycle extension parameter 


c x 


6.45 


5.0 


9.0 



regression (blue line), whose equation is given in the panel, and 
the set of expectation values for the model (y - x; red line). The 
correlation coefficient for the overall data set is approximately 
0.92 and the^ 2 per degree of freedom (i.e. reduced^ 2 ) is 0.16. 
Here, the reduced^ 2 value is computed for the worst case with 
the weight function (w) in Eq. (15) equal to 1. The slope of the 
regression (1.03 + 0.02) indicates a good reconstruction of the 
variability of the cycle amplitude, although the amplitude of the 
reconstruction is somewhat lower than the observations. In Fig. 
3b the distribution of the model error, defined here as the differ- 
ence between the model output and the observations is plotted, 
for the total magnetic flux. The mean error value is 0.55 x 10 14 
Wb and the standard deviation is about 9.03 x 10 14 Wb. The 
distribution is slightly asymmetric with a median value equal to 
1.42xl0 14 Wb. 

The scatter plot of the modeled versus the empirically recon- 
structed open flux is plotted in Fig. 3c. Annual averages are plot- 
ted. The regression (blue line) and the set of expectation values 
for the model (y = x; red line) are also displayed. The correla- 
tion coefficient is approximately 0.86 and the reduced x 2 , also 
computed with w = 1, is 0.31. The slope of the regression is 
0.86 + 0.05. For this parameter, the regression is biased by the 
high difference between the values observed and modeled for 
the maximum of cycle 19 and ascending phase of cycle 21. The 
Figure 3d shows the distribution of the model error for the open 
flux. The mean and median values for the error distribution are 
-0.17xl0 14 Wb and -0.19xl0 14 Wb, respectively. The standard 
deviation is 1.06 x 10 14 Wb. 

Since the open flux is distributed isotropically at 1 AU 
dLockwoo d et al. 2009a, and references therein), the radial field 
intensity at Earth is related to the open flux by 



Br(t) 



Anrl 



(18) 



where is the mean distance betw een the Earth and the Sun . 
The observed (red line; OMNI 2 data. lKing & Papitashvilil2005h 
and the modeled (blue line) radial component (\Br\) of the in- 
terplanetary magnetic field (IMF) are compared in Fig. |4^. The 
observations are averaged over Bartels rotations and cover the 
period from 1976 to 2009. FigureHJj presents the difference be- 
tween the observations and the modeled radial component of the 
IMF (blue line) and the standard deviation of values within a 
Bartels rotation (dotted green lines). We note that the observed 
values are systematically higher than the modeled values. The 
mean difference is about 0.25 nT with a standard deviation of 
0.68 nT (see Figs. 4d-e). In order to understand this discrep- 
ancy of 0.25 nT, we recall that to estimate the model parame- 
ter s we compared th e modeled output with the reconstruction 
by iLockwoodl (120091). which incorporate s a correction due to 
Kinematic effects dLockwood et al . 2009a b). This correction ef- 
fectively reduces the estimated open flux from interplanetary 



data. Consequently, the modeled radial component is system- 
atically lower than the interplanetary observations, that do not 
incorporate this correction. 

We test the null hypothesis that the mean value of |Z?r| on 
each Bartels rotation is the one computed by the model. Here, we 
apply the Student's test. Figure 4c presents the p-values within 
Bartels rotations, which are the probability of observing a value 
as extreme or more extreme of the test statistic value given by 



ttext — 



S 



(19) 



where (x) is the mean, \x is the modeled value, s is the sample 
standard deviation, and n is the sample size. For reference, the 
red dashed line displays the 5% significance level (a = 0.05). 
The red crosses in Fig. 4b indicate the periods at which we could 
reject the null hypothesis at the significance level of 5%. We note 
several periods in which the model values do not represent the 
averaged values over Bartels rotations. Long periods of discrep- 
ancy occur in the ascending phase of cycle 2 1 and the descending 
phases of cycles 21 and 23. 



3.2. Reconstruction of the solar magnetic fluxes since the 
Maunder Minimum 



The reconstruction of the total flux 



"total 



Oeph 



"open 



from 1700 to 2008 based on the group sunspot number is dis- 
played in Fig.|5^, while the reconstructions of the AR (blue line) 
and ER (green line) fluxes are plotted in Figure |5j5. Finally, the 
reconstruction of the open flux (green line) is shown in Fig.[5J:. 
For reference, the empirical reconstruction of the open flux (blue 
line) based on the geomagnetic aa-index is included in the plot. 
A clear secular trend in the ER and open flux leads to a secular 
trend in the total flux. 

As pointed out in Sect. 3.1.1, the modeled value of the total 
flux during the present minimum is approximately one third of 
the value observed during the previous minimum between cycles 
22-23. Furthermore, the modeled total flux during the present 
minimum is at the same level of the value returned by the model 
for the Dalton Minimum, a period of low solar activity lasting 
from approximately 1790 to 1830 (Figure 5a). The modeled to- 
tal flux during the present minimum is maintained entirely by 
the contribution of magnetic flux emerging in ERs belonging to 
cycle 23 since the flux emerging in ER for cycle 24 was not in- 
cluded due the lack of knowledge of parameters for this cycle. It 
must terefore be considered a lower limit. 

The group sunspot number, which extends from approxi- 
mately 1610 to the present, allows the solar magnetic flux to 
be reconstructed also prior to 1700. The model gives an open 
and total flux near zero through most of the Maunder Minimum 
(MM). This is a natural consequence of the coupling of the 
strength of the cycle of ER to that of sunspots, so that <p ep h turns 
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■ Modelled open flux 




140 




Fig. 2. (a) The open flux resulting from the present model (blue line) and the reconstruction based on the geoma gnetic aa-index from 
1904 to 2008 (black line). The dotted green line presents the open flux according to the model described by iKrivova et al.l ([2007 ). 
The square marks the value observed in 2008, the triangle, the modeled value. Cycle number is indicated at the bottom of the panel, 
(b) Rapidly (green line) and slowly (red line) evolving components of the open flux (blue line), (c) Fractional contribution of the 
rapidly (green line) and slowly (red line) evolving components to the total open flux. For reference, dashed red lines representing 
the value observed in 2008 are drawn in panels (a) and (b). 



out to be extremely weak in the Maunder Minimum. In addition, 
the length of the Maunder Minimum is longer than the decay 
time of the flux, so that in the model practically no flux survives 
from the pre-Maunder Minimum cycles until the end of the MM. 



4. Reconstruction of solar magnetic flux for the 
Holocene 

For studies aiming to isolate the Sun's influence on the Earth's 
climate time series of the solar magnetic and activity are needed 
as long as possible. Unfortunately, indices of solar activity such 
as the sunspot number, used here to reconstruct the AR, ER, to- 
tal and open magnetic flux, have been adequately recorded only 
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Fig. 3. (a) Scatter plot of the modeled versus the observed total flux, (b) Total Flux model error distribution, (c) Scatter plot of the 
modeled versus the reconstructed open flux. The cycle minima values (red squares) are indicated in the frame, (d) Open Flux model 
error distribution. 



since the invention of the telescope in the 17th century. They rep- 
resent the longest running time series of direct measurements of 
past solar variability. Cosmogenic isotopes provide estimates of 
solar activity that are less clean, in the sense that they are affected 
by the other quantities, such as the geomagnetic field and pro- 
cesses (e.g. climate variations and carbon cycle). However, these 
records e xtend to earlier times, cov ering periods up to thousands 
of years (Stuiver & Braziunas 1989|). 

By combining physics-based models for each of the pro- 
cesses connecting the isotope concentration in a relevant ter- 
restrial archive with solar activity, open flux and from it the 
suns pot number could be reco nstructed by Sola nki et al.l {2004) 
and lUsoskin et alj d2004l2007h . Since from cosmogenic isotope 
records the primary solar parameter that can be determined is the 
open flux, the problem faced here is opposite the one dealt with 
in Sect. 2 and 3. How to c ompute the SN from the open flux? 
This was first dealt with by Usos kin et alj d2002l) . In the recon- 
structions carried out so far, the solar open flux is linked with the 
sunspot number, AR and ER flux by inverting the model by de- 
scribing the evolution of the s olar surface magnet i c components 
for a given sunspot number dKrivova et al.ll200~7t ISolanki et alj 
2000, 2002). In this section, we update the previous estimate of 
the sunspot number, AR, and ER flux based on the inversion of 
the model described in Sect. 2.1. In addition, often only multi- 
year, e.g. decadal, averaged data are available. Hence we need to 
take this into account. 



4.1. Derivation of the magnetic flux model 

By differentiating Equation (O with respect to time, we obtain 



d(p 



dfopen dip, 



open 



dt dt dt 

Substituting Equations (HJl and (O in CO) , we find 



dip, 



open 



dt 



dt 



<ha_ 
^1 



<f>eph 
eph 



VQpen 



1 open 



' open 



After averaging over 10 years, we obtain 

d(popen\ i d( f> open 



dt 



dt 



'■'act 

IT - 



<Peph 
eph 



1 open 



1 open j 



(20) 



(21) 



.(22) 



Here, the symbol (...) denotes 10-year averaging. 

We assume that on decadal time scale the flux in active re- 
gions evolves in a steady state, i.e. the flux emerging in active 
regions is approximately equal to the decay due to several pro- 
cesses. In this way, we can write from Eq.Q] 



where 
1 1 



T act 



,-0 



(23) 



(24) 
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Fig. 4. (a) Modeled radial flux at 1 AU (blue line) in comparison with the measured radial interplanetary magnetic field component 
(OMNI data; red lines) averaged over Bartels rotations, (b) Difference between the measured and the modeled radial flux (blue line) 
at 1 AU. For reference, the 1-sigma value over Bartels rotations is also plotted (green lines). The read crosses (+) represent the 
values rejected in the hypothesis test (see text for the description of the test), (c) Student's test p-values over Bartels rotations. The 
red dashed line shows a = 0.05. (d) Scatter plot of the modeled versus observed radial flux at 1 AU. (e) Distribution of the difference 
between the measured and the modeled radial flux at 1 AU. 



The validity of this approximation can be tested by computing on a decadal scale is well founded. Substituting Eq. (0 into Eq. 
the ratio (<f> act ) I {e act ) from the model output, which should be ( [23) , we obtain 
close to T ac , ( 0.27 years). We found that ((f> acl ) / (e aa ) is approx- 
imately 0. 1 % higher than the value of r flc/ . Consequently, the as- 
sumption that the flux in active regions evolves in a steady state 

^nax,2 1 
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Fig. 5. Magnetic flux reconstruction from group sunspot number dHovt & Schattenlll998l) since 1700. (a) Total Magnetic Flux, (b) 
Active (blue) and Ephem eral (green) Region fluxes, (c) Modeled open flux (green line). The Open Flux reconstruction based on the 
geomagnetic index-aa by lLockwo od (2009) is plotted for reference (blue line). The dash red line represents the 2008 value of the 
Open Flux. 



Similarly, since 



eph 



14 hours is also very short compared 



to the cycle length, we obtain from Eqs. (f2]i and © 



(pephit) ~ T e ph£eph(t) ~ T eph X ^ ^™, X,l g l (f) 



where 
1 



Teph 



(26) 



(27) 



eph 



' eph 



The maximum emergence rate during the a ctivity cycle ; 
can be computed using Equation ([13). Following lUsoskin et al.l 
( 2007), we assume a linear relation between the amplitude of the 
solar cycle and the 10-year averaged sunspot number (R" wx '' = 



k < R „ >. The numerical value of k found by lUsoskin et al.l 
(120071) is 2.2 ± 0.4. The flux emergence in ER displays extended 
cycles, so that adjacent cycles partially overlap. For the set of pa- 
rameters presented in Table 1, adjacent cycles overlap in a way 
that a low 1 1 -year variability of the ER cycle is observed and the 
long-term trend is directly related to the 10-year averages of R g . 
At the maximum of an ephemeral region cycle, (f> ep h, we can set 
gi = 1, so that we obtain using Eq. (TT3T > 



/ max,i\ 
\£act ) 



max, 21 



max, 21 



f.21 



R" 



- k ( R z) 



(28) 
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by 



The ER flux averaged on a decadal time scale is then given 



max, 21 



(4>e P h)*T eph -^L-^kX(R g ) 



(29) 



We can also assume that the rapid open flux evolves in a 
steady state on a decadal time scale. In this case, 



open 



< <Kpen > 



(30) 



We found that the value of the ratio < (f> ac , >/< (f>' ope „ > is 
within approximately 0.5% of the value of ratio T r open /T 2 ct . 



Consequently, for a steady state, the term 



Eq. (|22) 



can be negleted. 

In order to compute the evolution of the open flux, we as- 
sume that 



fopen 



&(<f>open) 



At 



(31) 



\ dt I 
where 

(A<f>open) = {<f>open(tj + At)) - (<Popen{tj)) ■ (32) 

Substituting Eqs. (30) and (31) into d22l > we obtain 



(<fiopen) j +l \<t>apen). 



At 



1 / ^ 

1 1 open 

V act ' open 1 act J 



(<Peph) 

*>/+— T^,(33) 



eph 



where (4> pen) +1 = (<Po P en(tj + A?)}, {(popet). = (<f>open(t j)), and 



1 

At 



Replacing Eqs. (25) and (29) into (33) we obtain 

(<f>open) ((/>ope„) 

At Ti \ g/ J 

where, the constant c is given by 



1 



open 



> T 1 T s T 2 

V ' act ' open ' act 



Tact + 



T eph kX 



eph 



r*max,21 



g 



Rearranging, we get 

\ R g) ■ ~ a (<Popen) + b (tp pen). +i , 

where 
1 



and 



b = 



cAt 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



The 10-year averaged AR and ER can be retrieved by sub- 
stituting Eq. (37) in Eqs. (25) and (29), respectively. The rapidly 
evolving component of the open flux can then be obtained from 
Eq. (30). 



Empirical cumulative distribution function 




2 4 

Open Magnetic Flux (10 14 Wb) 

Fig. 6. Comparison of the empirical cumulative distribution 
functions, F(x), of the open flux. The blue line is the distri- 
bution of the modeled open flux based group sunspot number 
dHovt & Schatten|[l998l) . Red and black lines are distri butions of 
the open flux based on 14 C data dUsoskin e t al. 200% derived 



from paleo-geomagneti c reco nstructions of lYang et al.l ([2000) 
and iKorte & Constable! d2005l) . respectively. The 2008 annual 
average is indicated by the red circle. 



4.2. Reconstructed magnetic flux through the Holocene 

Our reconstruction of past solar activity relies on the decadal es- 
timate of the so lar open flux from m easu rements of 14 C as wa s 
done earlier by ISolanki et all d2004 and lUsoskin et al.l d2007l) . 
As discussed by lUsoskin et alJ d2007l) . the estimate of the open 
flux depends on the k nowledge of the tem poral evolution of 
the geomagnetic field. Usos kin et alJ (I2007I) presen ted two re- 
constr uctions b ased on the paleomagne tic models by Yan g et all 
(120001) and IKorte & Constable! <|2005). The first one extends 
through the whole Holocene while the second one reaches back 
around 7000 years. In Fig. [6] a comparison between the empir- 
ical cumulative distribution functions of these two reconstruc- 
tions based on I4 C data and the one based on the telescopic 
sunspot record is presented. In order to compare the distribu- 
tions of the open flux based on 14 C, we have employed just the 
period over which the time series overlap. A s the geomagnetic 
dipole moment of Korte & Constable ( 2005b is systematically 
lower than that obtained by Yang et al] d2000l) . a systematically 
higher open flux is obtained. We note that the distribution of 
the open flux estimated from sunspot number is closer to the 
reconstruction using the Korte & Constable (2005) geomagnetic 
dipole momen t. It suggests that if the p aleo-geomagnetic recon- 
structions of IKorte & Constable ( 2005) is close to the real evo- 
lution of the magnetic field, the values of the open flux since 
Maunder Minimum are not unusual comparing to the values ob- 
served during the Holocene. I f, howe ver, the paleo-geomagnetic 
reconstruction by Yang et alJ (l2000l) is closer to the real evo- 
lution of the magnetic field, the values of the open flux since 
the Maunder Minimum are unusually high compared to the 
Holocene. Furthermore, we note that about 25% of the modeled 
values of the open flux from these two reconstructions are be- 
low the value observed in the present minimum based on the 
IKorte & Constable! (120051) reconstr uction, while abo ut 45% of 
the values of the open flux based on lYang et al.l ([2000) are below 
the present value. Note that the value for the present minimum 
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Fig. 7. Lon g-term total magnetic flux reconstruction from 14C data. The blue line is the reco nstruction based on the open flux 
obtained by Usosk irTet al.l (|2007)using paleo-geomagnetic data from iKorte & Constable! ((2005). The green line is the reconstruc- 
tion of the total magnetic flux (10-year running means) based on the group sunspot number since 1700 AD. Panel (a) shows the 
reconstruction from 500 AD to present while panel (b) shows the reconstruction since 5000 BC. 



is a yearly value, while the curves based on the reconstructions 
from 14 C are decadal averages. 

A comparison between the reconstructions of the sunspot 
num ber obtained using Equation (37) and the one obtained 
by lUsoskin et al.l d2007l) reveals a good correspondence (R « 
0.96). We note that the reconstructions obtained based on the 
two approach e s are quite similar, with the reconstruction from 
lUsoskin et al.l (120071) having slightly lower values during high 
activi ty. Both reconstruct ions are bas ed on the open flux e s timate 
by lUsoskin et all J2007) using the IKorte & Constable! d2005l) 
model. The main difference between the two models is the rela- 
tionship between the open flux and sunspot number. In this work, 
we distinguished the fast and slowly evolvin g component s of the 
open flux, while in the reconstruction by Usos kinet alj d2007l) 
such a separation is not made. 

We obtain the reconstruction of the total magnetic flux from 
the estimate of the magnetic flux in A R (Eq. (25)), ER (Eg . 
(29)), and the open flux provided by Usosk inet al.l (|2007). 
Figure [7] presents the o btained reconstruction based on the 
Kort e & Constable] (120051) geomagnetic model. In Fig. 7a, we 
concentrate on the reconstruction since 500 AD (blue line). For 
reference, the total flux estimated from sunspot data is presented. 
The reconstruction since 5000 BC is shown in Figure 7b. The to- 
tal flux is required to compute the irradiance, which will be the 
topic of a forthcoming paper. 



5. Concluding remarks 

In the present p aper, we have cons idere d an extension of the sim- 
ple model of ISolanki etail d2002l) and lKrivova et al.l (|2007|) de- 
scribing the evolution of the Sun's open and total magnetic flux. 
We have shown that by considering separately a rapidly evolv- 
ing and a slowly evolving component of open flux we obtain a 
greatly improved agree ment with of the s olar open flux recon- 
structed since 1904 by iLockwoodl d2009l) and a reasonable al- 
though not perfect description of the OMNI data as well as. The 
main improvement provided by this version of the model is in the 
reproduction of the cyclic variation of the open flux, including 
the amplitudes of individual cycles. 

The rapidly decaying open flux is most likely harbored in 
small coronal holes associated with ARs or decaying ARs, while 
the slowly decaying open flux is associated with the polar coro- 
nal holes, but also with open flux at low latitudes. 

We found that approximately 25% of the modeled open 
flux values since the end of the Maunder Minimum are lower 
than the low observed during the present minimum (i.e. in 
2008), which suggests that the present solar minimum condi- 
tions are not exceptional in terms of the heliospheric magnetic 
flux. We noted that the same amount is observed in the recon- 
struction of the open flux by lUsoskirTe t al. (2007) based on the 
IKorte & Colrsta ble (2005]) model, while about 45% of the open 
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flux values are lower than the valu e during the curre nt minimum 
in the reconstruction based on the Yang et aT] d2000h model. 
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